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Abstract. During the gravitational core collapse of a massive progenitor star which may give rise to at least a class of 
gamma-ray bursts (GRBs) associated with supernovae, a stellar core rapidly passes through a short yet important phase of 
neutronization, producing a huge amount of energetic neutrinos and photons which contribute to the total pressure within 
the progenitor core. The collection of neutrinos, photons and gas materials together may be approximated as a fluid with a 
polytropic index y = 4/3 under the action of self-gravity. With a substantial generalization and using analytical and numerical 
methods (Lou & Cao 2008), we recently constructed and examined various self-similar solutions to describe collapses, 
rebound shocks and flows systematically in a y= 4/3 polytropic gas mixture with spherical symmetry, and compare our results 
with those of Goldreich & Weber (1980). It is also possible to construct central void solutions without or with shocks. Various 
features and characteristics of this nonlinear relativistically hot gas dynamics, including asymptotic and exact solutions, are 
presented. This more general polytropic model analysis provides the dynamic basis of understanding the evolution of rebound 
shocks in supernovae (SNs) and the results may be also utilized to benchmark hydrodynamic simulations. 
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INTRODUCTION 

Long gamma-ray bursts (GRBs) are generally thought to originate from demises and explosions of massive stars. They 
release the gravitational energy by drastic core collapses and shock outbursts, producing characteristic radiations in 
different energy bands, including gamma rays, X-rays and energetic neutrinos. While the key mechanism of GRBs 
still remains mysterious, various observations offer valuable clues to these extremely explosive "cosmic fireworks". 
Especially since the discovery of the precise match in positions between GRB980425 and SN1998bw, several GRBs 
are clearly associated with at least one class of supernovae (Type Ibc) which are related to deaths of old massive stars. 
Three other such associations GRB030329 and SN2003dh, GRB031203 and SN20031w, GRB060218 and SN2006aj 
together with a recent X-Ray Outburst (XRO)-SN connection between XRO 080109 and SN 2008D (e.g., Woosley & 
Bloom 2006; Soderberg 2008) were reported within the last decade. For long GRBs, dynamics of stellar core collapses 
and rebound shock evolution in the context of SN explosions are important to both theories and observations of GRBs. 

With certain assumptions about the state of collapsing stellar cores, hydrodynamics and magnetohydrodynamics 
(MHD) may be adopted to model these stellar collapses with or without shocks (e.g., Lou & Wang 2006, 2007). Based 
on gas properties, a self-similar dynamics may emerge in these astrophysical processes. First, to simplify the problem, 
a quasi-spherical symmetry is presumed and only self-gravity and thermal pressure are taken into account in this semi- 
analytic model of a relativistically hot polytropic gas. A simple yet useful assumption is the conventional polytropic 
state function, i.e., P = KpY where P and p are pressure and density, /is the polytropic index and K is a global constant. 
Based on such an assumption, valuable information about the stellar structure, solar physics, star formation and SN 
explosions may be obtained (e.g., Chandrasekhar 1939; Shu 1977; Lou & Shen 2004; Suto & Silk 1988; Lou & Wang 
2006, 2007). A simple argument on total energy of a static polytropic sphere demonstrates that y > 4/3 cases are stable 
while 7 < 4/3 ones are unstable (Chandrasekhar 1939). Thus most researchers focus on 1 < y < 4/3 cases. 

We emphasize that the critical case of y = 4/3 is important in astrophysical contexts. Statistical physics illustrates 
that in a temperature much lower than the Fermi energy, the state of particles whose rest mass is negligible is a 
y = 4/3 polytrope. Meanwhile, numerical simulations show that y = 4/3 polytrope is a sensible approximation for 
the central state of a compact object. In addition, the pressures of neutrinos and of photons, both crucial in collapsing 
stellar cores, may also be approximated by a y = 4/3 polytropic state. It turns out that the conventional polytropic 
assumption limits y = 4/3 situation severely in a self-similar model. In a classic paper on y = 4/3, Goldreich & 
Weber (1980; GW hereafter) considered homologous collapse and concluded that when the pressure decreases by a 
fraction of no more than 2.9%, a homologous core collapse would occur in the stellar interior, which is considerably 



lower than a value of 26% from numerical simulations of Bethe et al. (1979). GW tried to reduce this difference by 
introducing an inner core in a progenitor. Yahil (1983) performed a polytropic analysis and treated GW results as a 
limit of y-> (4/3)~. 

In contrast to the conventional gas of a constant K, a more general polytropic gas still obeys 

for the specific entropy conservation along streamlines with a variable P/p? (e.g., Cheng 1977, 1978). Clearly, the 
conventional polytrope is only a special case. Nonetheless, mathematical derivations become more involved if the 
more general polytrope is adopted in the model analysis. 

Recently, Lou & Cao (2008) used the self-similar transformation (Fatuzzo et al. 2004) for a more general polytrope 
of 7 = 4/3 (see eqn (fl]i). Based on different values of a scaling index a, the problem is divided into three classes. Some 
essential procedures, analyses, and results are summarized in the following. 



MODEL FORMULATION AND MAIN RESULTS 

In spherical polar coordinates (r, 0, (j)) with spherical symmetry, physical variables are functions of radius r and time t. 
Conservations of mass and momentum together with the local specific entropy conservation along streamlines govern 
the dynamic evolution of a relativistically hot gas of y = 4/3. Without shocks, the problem is invariant under the time 
reversal transformation. Thus, an expansion solution can also describe a collapse solution by such a transformation. 
In order to solve hydrodynamic partial differential equations (PDEs), a self-similar transformation is introduced to 
reduce PDEs to ordinary differential equations (ODEs), namely 
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where p,M, u and P are respectively density, enclosed mass, radial velocity and pressure, a is a scaling index in such a 
transformation, G is the gravitational constant and A is a dimensional parameter to make the independent self-similar 
variable x dimensionless. Then the governing hydrodynamic PDEs are reduced to nonlinear ODEs 
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For considerations of dynamic processes in the core of a massive star and especially rebound shocks of SNe associated 
with GRBs, we set y = 4/3 to approximate the state of gas materials in the compact central core. Generally speaking, 
we require a being less than zero. Importantly, equation (0 implies a division of all cases into three classes, depending 
on the value of a being greater than, equal to or less than —2/3. In the following analysis, we deal with these classes 
separately and obtain respective results. Comparing to GW, our analysis to the a = —2/3 class represents a substantial 
generalization. 

The case of a = —2/3. For a = —2/3, the reduced specific entropy conservation (O becomes automatically 
satisfied. Meanwhile from equations ([3])— ©, one can readily derive the following two relations 
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It is taken that the pressure P is proportional to p 4 / 3 , and we need to know the proportional coefficient as a function 
of (r, t) to complete the problem. Physically, log(P/p 4 / 3 ) is proportional to specific entropy s(r, t) in a polytropic 



gas. Once this distribution s(r, t) is known, the gas dynamics is determined. A more general case in self-similar 
transformation is to allow p = g(x)a 4 ^ where g(x) is a function of x to be specified for a self-similar flow. For a given 
g(x), equation (jT) can be solved with proper initial or boundary conditions. In general, velocity cannot diverge at large 
radii. Therefore a boundary condition to this problem is that an outer boundary of zero density should exist. 

The first cut is to adopt a constant specific entropy s(r, t) everywhere, i.e., g(x) = 1, as GW did. We confirmed the 
main results of GW. An essential result is that the central density has a critical minimum below which the solution 
of equation (O has no vanishing a at a finite radius. This critical value of a is ~ 101.88, consistent with GW work. 
In the context of a stellar collapse, GW result is somewhat limited in certain aspects. Take an old star in a relatively 
stationary state as described by the Lane-Emden equation (e.g., Chandrasekhar 1939). As the central nuclear burning 
becomes insufficient, the whole system rapidly evolves into a homologous collapsing or pre-collapsing phase. GW 
estimated the largest pressure reduction fraction (denoted by r c here) for a homologous collapse (prior to a rebound) 
and obtained a value of ~ 2.9% under their assumption for a polytropic sphere of y — 4/3. In contrast, numerical 
simulations of Bethe et al. (1979) show that this value r c in supernovae could be substantially larger, about 26%. To 
reconcile, GW proposed that there is an inner and smaller core which obeys their results as the pre-collapsing phase. 

To generalize GW analysis for a homologous collapse, we allow a fairly arbitrary form of g(x) in our model and thus 
accommodate a broad class of solutions for the density profile. Such a variable s(r, t) may lead to a consistent result 
with numerical simulations (e.g., Bethe et al. 1979). In our numerical exploration, we choose g(x) = 1/(1 + ex) with 
£ being a coefficient to measure the range of g(x) variation. According to our experiments, when £ = 0.1, the pressure 
reduction r c can reach ~ 10% and when £ = 0.3, r c is ~ 26% as given by Bethe et al. (1979). Given idealizations, these 
explorations do reveal a simple fact that a variable distribution of specific entropy s(r, t) can be extremely important 
and we could model this effect in a self-similar manner. 

The case of a < —2/3. In these cases, equations (0) and © lead to a new form of state function below 

P = c m^, g = g|±£=rt (8) 

(3a + 2) 

where Co is a constant. Most importantly, the variable specific entropy \og(P/pY) is connected with the enclosed mass 
mix) in a power-law form. The conservation of mass guarantees the conservation of specific entropy. For q = and 
thus a = 7 — 2, we have a conventional poly trope; if we further require 7=1 and thus a = — 1, we then have an 
isothermal gas (Shu 1977; Lou & Shen 2004). For 7 = 4/3 under consideration, q = 2/3, i.e., p = Co/« 2//3 0! 4//3 , which 
is a completely different family of problem in contrast to P = Jfp 4 / 3 . Note that for 7 =/= 4/3 cases, the constant Co 
can be absorbed by readjusting the A parameter in the transformation but for 7 = 4/3, constant Co is independent and 
characterizes the local sound speed. 

It is easy to derive two coupled nonlinear ODEs for v(x) and a(x) from equations dill—©. Through analysis, we 
obtain all kinds of relevant solutions both asymptotically and numerically, with two globally analytical solutions. 
Numerical simulations are needed to decide the paths of various self-similar evolution. 

In our analysis, we compare with known results, finding more general counterparts of known situations, such as 
the central free-fall solution (Larson 1969; Penston 1969), the expansion wave collapse solution (EWCS; Shu 1977), 
envelope expansion with core collapse solutions (EECC; Lou & Shen 2004) and quasi-static solutions (Lou & Wang 
2006, 2007). Furthermore, two new types of asymptotic solutions are also found. Corresponding physical meanings 
are discussed for possible processes in astrophysics. To enlarge the range of astrophysical applications, shock waves 
are also included. For shock conditions, mass, momentum and energy should be all conserved across a shock front. 
Numerical solutions connecting different types of asymptotic solutions are constructed in our model, with relevant 
physical scenarios outlined (see Lou & Cao 2008 for details). 

Most interestingly, we can construct central void solutions without or with shocks in the framework of hydrody- 
namics involving thermal pressure and self-gravity. In a background of Einstein de Sitter expanding universe, Fillmore 
& Goldreich (1984a, b) considered a collection of collisionless particles with self-gravity. In the presence of certain 
kinds of perturbations, voids can be generated surrounding the center of such a flow of collisionless particles. Simi- 
larly in our model, although without an expanding background, thermal pressure also has a tendency of expansion and 
thus provides a driving mechanism to produce expanding voids under favorable conditions. Observationally, voids are 
ubiquitous in various astrophysical contexts, such as supernova remnants, molecular clouds, hot bubbles, superbubbles, 
planetary nebulae and galaxy clusters (see Hu & Lou 2008 and Lou & Hu 2008 for more details). 



The case of a > —2/3. Since the enclosed mass cannot be negative, a > —2/3 requires v > —ax in equation (0). 
However, as we cannot obtain any solution with an outer boundary of finite radius, velocity should converge as x 



goes to infinity. However, these two requirements are incompatible with each other. The requirement of v > -ax will 
certainly lead to a divergent velocity at large x. These cases are thus excluded in our model consideration. 

SUMMARY AND CONCLUSIONS 

Both analytic and numerical self-similar solutions have been explored for a more general polytropic gas with 7 = 4/3 
(Lou & Cao 2008; Wang & Lou 2008). We have directly extended the classical analysis of GW for a conventional 
polytropic gas of spherical symmetry with K = 1, giving a plausible resolution to the discrepancy between theoretical 
calculations (GW) and numerical simulations (Bethe et al. 1979). We then study distributions of specific entropy in a 
certain system or infer these from other available information on the basis of our theoretical knowledge. Once such an 
entropy distribution is known, the self-similar dynamical evolution of a homologous core collapse can be calculated 
readily from our model. Another family of solutions are also discussed with a new state function of P « M 2 / 3 p 4 / 3 . 
In addition to counterparts of various previously know types of analytic and asymptotic solutions at y ^ 4/3, we find 
some new ones and discuss their physical implications. In Lou & Cao (2008), we attempt to analyze the self-similar 
dynamics of a y = 4/3 polytropic gas systematically and hope that it would provide valuable and useful clues to the 
study of compact objects and high-energy explosive phenomena, such as shock breakouts of SN explosions and thus 
GRBs. 

There are several physical processes, such as nuclear reactions, radiation pressure, neutrino transportation, general 
relativity, rotation and magnetic field effects, that are not included in our preliminary model analysis. However, given 
approximations and idealizations of our model, we hope that the results of this theoretical hydrodynamic analysis catch 
certain essential features of large-scale shocks in a relativistically hot gas. 
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